Sequential melting of charmonium states in an expanding Quark Gluon Plasma and 
J/ if) suppression at RHIC and LHC energy collisions 
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We have developed a hydrodynamic model to study sequential melting of charmonium states in 
an expanding QGP medium. According to the initial fluid temperature profile, J/^'s are randomly 
distributed in the transverse plane. As the fluid evolve in time, the free streaming J/ip's are 
suppressed if the local fluid temperature exceeds a critical temperature. PHENIX data on the 
centrality dependence of J/tp suppression in Au+Au collisions at mid-rapidity are explained by 
sequential melting of the charmonium states, Xc, and J/ip, in the expanding medium. The 



critical temperatures T 



j/i, 



2.09Tc and 



T 



1p/ 



l.lTc agree with lattice motivated 



calculations. The feed-down fraction F depend on whether the cold nuclear matter effect is included 
or not. It changes from F = 0.3 with cold nuclear matter effect included to F = 0.5 when the effect 
is neglected. Model fails to reproduce the PHENIX data on the centrality dependence of J/^ 
suppression in Cu-I-Cu collisions at mid-rapidity, indicating that the mechanism of J/tp suppression 
is different in Au-|-Au and in Cu-|-Cu collisions. We also use the model to predict for the centrality 
dependence of J/ip suppression in Pb-|-Pb collisions at LHC energy, -^8=5500 GeV. In LHC energy, 
J/i/i's are more suppressed in mid central collisions than in Au+Au collisions at RHIC energy. 



PACS numbers: PACS numbers: 25.75.-q, 25.75.Dv 



I. INTRODUCTION 

In relativistic heavy ion collisions J/^ suppression has 
been recognized as an important tool to identify the pos- 
sible phase transition to quark-gluon plasma. Because 
of the large mass of the charm quarks, cc pairs are pro- 
duced on a short time scale. Their tight binding also 
makes them 'nearly' immune to final state interactions. 
Their evolution probes the state of matter in the early 
stage of the collisions. Matsui and Satz [l|, predicted 
that in presence of quark-gluon plasma (QGP), binding 
of a cc pair into a J/ip meson will be hindered, leading 
to the so called J/ip suppression in heavy ion collisions 
[H . Over the years, several groups have measured the 
J/^ yield in heavy ion collisions (for a review of the data 
prior to RHIC energy collisions, and the interpretations 
see Refs. 0, Q)- In brief, experimental data do show 
suppression. However, this could be attributed to the 
conventional nuclear absorption, also present in pA col- 
lisions. 

PHENIX collaboration has undertaken the task to 
characterise medium effect on J/ip production in nuclear 
collisions. They have measured J/tp yield in p+p colli- 
sions at RHIC and obtained the reference for basic in- 
variant yield [1, H, 0. Measurements of Jf-ip production 
in d-|-Au collisions [a, 0| give reference for cold nuclear 
matter effects. J/tp production in d-|-Au coUisions are 
consistent with cold nuclear matter effect quantified in a 
Glauber model of nuclear absorption with (Tabs = 2 ± 1 
mb [8]. Cold and hot nuclear matter effects are stud- 
ied in Au-|-Au and Cu+Cu collisions, where yields are 



measured as a function of collision centrality [T^, [HI ■ 
In Au-I-Au and Cu-|-Cu [ll[ collisions, data are 

taken at mid-rapidity (|y| < .35) and at forward rapidity 
(1.2 < y < 2.2). J/ip's are more suppressed at forward 
rapidity than at mid rapidity. It was also noted that 
at comparable participant number, J/t/^'s are suppressed 
similarly in Au+Au and in Cu+Cu collisions 

At RHIC energy, it has been argued that rather than 
suppression, charmonium's will be enhanced (l^ . [isj . 
Due to large initial energy, large number of cc pairs will 
be produced in initial hard scatterings. Recombination 
of cc can occur enhancing the charmonium production. 
Both the PHENIX data on J/ip production in Au+Au 
and in Cu+Cu collisions, are not consistent with mod- 
els which predict J/ip enhancement [H, The cold 
nuclear matter effect, quantified by the Glauber model 
of nuclear absorption with J/ V'-nucleon absorption cross- 
section <7abs = 2±1 mb, is consistent only with peripheral 
Cu+Cu collisions. In all the centrality ranges of Au+Au 
collision, or in central Cu+Cu collisions, suppression is 
beyond the cold nuclear matter effect. 

Blaizot et aZ.fl^. [TH proposed a phenomenological 
model, called the threshold model, to describe J/ip sup- 
pression in Pb+Pb collisions at SPS energy. In the 
threshold model, to mimic the onset of deconfining phase 
transition above a critical energy density and subse- 
quent melting of J/ip, J/ip suppression is linked with 
the local transverse density. If the local transverse den- 
sity at the point where J/ip is formed exceed a criti- 
cal or threshold value, J/ip's are melted. Recently, in 
the threshold model, we have analysed the PHENIX 
data on the centrality dependence of J/ip suppression in 
Au+Au/Cu+Cu coUisions at mid-rapidity 16, 17, isj. 
Sequential melting of charmonium states, Xci i^' and 
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J/ip above threshold density. 
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explains the PHENIX data on J/-0 suppression in mid- 
rapidity Au-I- Au collisions [3 ■ It was also observed that 
the feed-down fraction F, from higher states Xc and ip', 
depends on the quantum of nuclear absorption. Equiv- 
alent fit to the data could be obtained by increasing 
the fraction F and decreasing the J/^/'-nucleon absorp- 
tion cross-section or the vice- versa. The threshold model 
ignore the expansion of the medium. Also, in the thresh- 
old model, while it is assumed that the threshold density 
is proportional to the critical energy density above with 
charmonium states melt, the exact relation between them 
is rather obscure. 

Recently Gunji et al.fl^ analysed the PHENIX mid- 
rapidity data on the centrality dependence of J/iJj sup- 
pression in Au+Au collisions . They developed a 
"Hydro+J/^" model. The QGP fluid evolves in 3+1 di- 
mensions. At the initial time, in accordance to the fluid 
temperature, J/ip^s are randomly distributed in the fluid. 
As the fluid evolve, free streaming J/ip^s are melted if 
the local fluid temperature exceeds a critical value. The 
experimental J/^ suppression pattern in mid-rapidity 
Au-|-Au collisions is well explained by sequential melting 
of Xc ip^ and J/^ in the dynamically expanding fluid. 
The estimated melting temperatures, Tj/^ — 2.02Tc, 
T^^ — T^, — 1.22Tc are in agreement with the lattice 
motivated calculations [13]. The fraction of the higher 
states (xc + ■0') is estimated to be F = 0.3. It may be 
mentioned that experimentally, the feed-down fraction F 
is largely uncertain [2l| . Measurements are available over 
a wide range of energy -^5=8. 5-1800 GeV. Measured val- 
ues show considerable variation, F — 0.15 — 0.74. The 
estimated fraction _F=0.3 is well within the largely uncer- 
tain range of measurements. Incidentally, measurement 
at an energy comparable to RHIC energy y^=200 GcV 
is not available. 

In the present paper, we have developed a hydrody- 
namic model for J/^ suppression in heavy ion collisions. 
The aim was to verify the results obtained by Gunji et 
aLp^ and to use the model to predict for the suppres- 
sion pattern at Pb+Pb collisions at LHC energy. The 
model is similar to that of Gunji et al.fl^, though there 
are some differences in details. Gunji et al. 19] solved 
hydrodynamic equations in 3+1 dimensions. The ini- 
tial fluid energy density in the transverse plane was pa- 
rameterised in proportion to the Glauber model calcu- 
lation of hard collisions. We have solved the hydrody- 
namic equations in a 2-1-1 dimensions assuming longi- 
tudinal boost-invariance. The model is thus limited to 
J/V' production at mid-rapidity only. The initial fluid 
energy density in the transverse plane is assumed to be 
proportional to the Glauber model calculation with 75% 
soft collisions and 25% hard collisions [13]. As it will 
be shown later, PHENIX data [9., j3] on the centrality 
dependence of J/'0 suppression in Au+Au collisions arc 
reproduced in the model with melting temperatures and 
feed-down fraction close to the values obtained by Gunji 
et al.^^. We have also analysed the recently published 
Cu+Cu data [ll|. Hydro+J/V-' model with melting tem- 



peratures extracted from the analysis of Au+Au data, 
is not consistent with the J/ip suppression in Cu+Cu 
collisions. The model produces more suppression than 
required by the data. Clearly, J/?/; suppression mecha- 
nism in Au+Au and Cu+Cu collisions is not identical. 
We also use the model to predict for the centrality de- 
pendence of J /'ip suppression in Pb+Pb collisions at LHC 
energy (y^=5500 GeV). In LHC energy collisions, J/^'s 
are more suppressed than in Au+Au collisions at RHIC. 

The plan of the paper is as follows: in section II, we 
briefly describe the hydrodynamic model. The explicit 
mechanism followed to suppress J/ip^s in the expanding 
fluid is also discussed in section II. In section HI, we have 
analysed PHENIX data on the centrality dependence of 
J/V' suppression in Au+Au collisions and obtain the pa- 
rameters of the model. We also analyse the PHENIX 
data on the centrality dependence of J/^p suppression in 
Cu+Cu collisions. In section IV we give predictions for 
J/tp suppression in Pb+Pb collisions at LHC energy col- 
lisions. Lastly, summary and conclusions are drawn in 
section V. 



II. THE HYDRODYNAMIC MODEL FOR J/ip 
SUPPRESSION 

A. Hydrodynamic model for QGP evolution 

Details of the hydrodynamic model used here can be 
found in [2^ ]. In (2^ . Kolb and Heinz, assuming lon- 
gitudinal boost-invariance, solved the energy-momentum 
conservation equation cJ^T'"' = 0, in 2+1 dimensions. 
Hydrodynamic models require energy density, fluid ve- 
locity distributions at the initial time r^. In the 
initial energy density of the fluid in the transverse plane 
was parameterised as, 

£(b, X, y) = eo[0.75Npartih, x, y) + .25iVco;/(b, x, y)] (1) 

where A^part(b, x, ?/) and iVco/;(b, a;, y) are the transverse 
profile for the participant number and binary collisions 
number in an impact parameter b Au+Au collision. The 
initial fluid velocity was assumed to be zero, i'a;(b, x, y) = 
Vy{h,x,y) = 0. The constant depend only on the col- 
lision energy and not on centrality of the collisions. The 
initial time and the constant chosen to reproduce the 
Pt distribution of identified particles in central Au+Au 
collisions. For b=0 Au+Au collisions, it correspond to 
central energy density, £=30 GeV/ fm? or central entropy 
density S'i„i=110/m~'^, at the initial time Ti=0.6 fm/c. 
Kolb and Heinz [22, used an equation state (EOS-Q) in- 
corporating 1st order phase transition with critical tem- 
perature rc=164 MeV. The quark phase was modeled by 
the bag equation of state for u,d,s quarks and gluons. 
For the hadronic phase, resonance hadron gas equation 
of state was used. The bag constant was obtained by 
using the Maxwell construct at the critical temperature 
Tc=164 MeV. 



3 



2.5 



(a)T(x,y=0), b=2.3 fm (b)T(x,y=0), b=7.4 fm 




Au+Au 

Tc=0.164 GeV 



1.0 
2.5 - 

2.0 - 



-5 5 

I 

(c)T(x=0,y), b=2.3 fm 



(d)T(x=0,y), b=7A fm 




FIG. 1: In panels (a) and (c), the evolution of local tem- 
perature T{x,y = 0) as a function of x and T{x = 0,y) as 
a function of y, in b=2.3 fm Au+Au collisions are shown. 
In the panels (b) and (d) the same are shown for b=7.4 fm 
Au+Au collisions. In each panel, the different lines, from top 
to bottom correspond to ri=0.6,1.6,2.6 and 3.6 fm. 



Evolution of QGP fluid in the model is described in 
detail in [l^]- For completeness, in Figllfa) and (b), we 
have shown the evolution of local temperature T{x, y = 
0), in unit of the critical temperature Tc, in b=2.38 fm 
(panel (a)) and b=7.4 fm (panel (b)) Au+Au coUisions. 
In FiglUc) and (d) local temperature T{x = 0,y) as a 
function of y is shown. The different lines (from top 
to bottom), in each panel, corresponds to time Ti— 0.6, 
1.6, 2.6 and 3.6 fm/c. Lattice motivated calculations 
indicate the just at the critical temperature Tc, all the 
charmonium states are not dissolved [13] ■ The ground 
state J /ipilS) can survive in QGP environment up to a 
temperature Tj/^ « ^.ITc- The excited states Xci^P) 
and tjj'{2S) on the other hand cannot survive hot QGP. 
They are dissolved in much cooler QGP, T^^ « 1.2Tc 
and T^i « l.lTc. From FiglU it is obvious that while 
the excited states Xc and ij^f will be melted in b=2.3 and 
7.4 fm collisions, the ground state can be melted only in 
b=2.3 fm Au+Au collisions. The peak temperature in a 
b=7.4 fm collision is well below the melting temperature 
of the ground state J/i^- Indeed, without any detailed 
calculations, from the initial peak temperatures in dif- 
ferent centrality collisions, one can very well bound the 
melting temperature for J/^ and the states Xc and ■(/;'. 
Initial temperature T{x,y = 0) of the fluid in b=: 2.3, 
4.1, 5.2, 6.2, 7.4, 8.7, 10.5 and 13.0 fm Au+Au collisions 
are shown in FigO Roughly they corresponds to 0-5%, 
5-10%, 10-15%, 15-20%, 20-30%, 30-40%, 40-60% and 60- 
93% Au+Au collisions. In 0-5% centrality collision, peak 
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FIG. 2: Initial temperature, in unit of the critical tempera- 
ture, of the QGP fluid in Au+Au collisions at impact param- 
eter b= 2.3, 4.1, 5.2, 6.2, 7.4, 8.7, 10.5 and 13.0 fm (from top 
to bottom). They corresponds to 0-5%, 5-10%, 10-15%, 15- 
20%, 20-30%, 30-40%, 40-60% and 60-93% centrality Au+Au 
collisions. 



temperature is only 2.16rc. As the collision centrality 
decreases, peak temperature decreases and in the most 
peripheral (60-93% centrality) collisions peak tempera- 
ture is only 1.16Tc. As shown in Figdl in all these cen- 
trality ranges of Au+Au collisions J/tp^s are suppressed 
Presumably, the suppression in 60-93% centrality 
collisions is due to melting of the higher states Xc and ip' 
only. One then obtain an upper bound for the melting 
temperature of the states Xc and ip/, T^^,T^„ < 1.16Tc. 
Similarly, in 0-5% centrality collisions, if the ground state 
J/'f/^'s are dissolved then melting temperature can be 
bounded from above, Tj/^ < 2.16Tc. The bounds on 
the melting temperature are close to the lattice results 
[20t or to the values obtained by Gunji et al.^^ from the 
analysis of PHENIX data. 



B. J/i/> suppression in the QGP fluid 

To obtain the survival probability of J/^/^'s in an ex- 
panding QGP fluid, we proceed as follows: at the initial 
time Ti — 0.6 fm/c, according to the initial spatial distri- 
bution of the fluid temperature, we randomly distribute 
J/^j in the transverse plane x — y. The transverse mo- 
mentum of initial J /ip are distributed according to the 
power law A/{1 + {pt/B)^)^, which rather well describe 
the invariant distribution of measured J/ip's in p+p col- 
lisions The initial J/V''s are assumed to follow a 
free streaming path unless dissolved in the medium. To 
follow the path, to each J/ip, we assigned a random direc- 
tion vector. It may be mentioned that while assumption 
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FIG. 3: (color online) The black filled circles are PHENIX 
data on the centrality dependence of J/tp suppression in 
Au+Au collisions. The black line is the cold nuclear 
matter effect, computed in a Glauber model with J/tp- 
nucleon absorption cross-section o-abs=2 mb. The dotted lines 
with filled colored symbols are J/tp survival probability in 
the present "Hydro+ J/?/)" model, for melting temperature 
Tj/^/rc=l. 2, 1.4, 1.6, 1.8,2.0 and 2.2 respectively 



of free streaming do not affect the centrality dependence 
of J/^ suppression , it will affect the px distribution . 
Experimentally, px distribution of J/ijj^s are broadened. 
The free streaming assumption will not reproduce the px 
broadening of J/t/j's. 

Following Gunji et oZ.[19i]. the survival probability of a 
J/4> inside the expanding QGP is defined as, 



Sj/^{t) = exp 



r<i,,(T(x(T')))dT' 



(2) 



where T(x) is the temperature of the fluid at the trans- 
verse position X, Tdis{T) is the decay width of J/^ at 
temperature T. ti is the initial time for hydrodynamic 
evolution. We continue the evolution till the freeze-out 
temperature Tf=130 MeV. For the decay width of J/V'j 
we made the simple choice, 
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FIG. 4: (color online) The black filled circles are PHENIX 
data on the centrality dependence of J/tp suppression in 
Au+Au collisions. The black and red lines are best fit to 
the data in the model, with parameter values indicated in the 
figure. 



B-like states can exist in sQGP. They can provide reso- 
nant cross section for heavy quarks. Recently van Hees et 
aL[25|, in a relativistic Langevin approach showed that 
resonant interactions play an essential part in thermalisa- 
tion and collective flow of charm and bottom quarks. D- 
and B-like states in sQGP can also have direct impact on 
J/f/' production, facilitating regeneration. Survival prob- 
ability (Eq|21) also neglects regeneration of charmonium 
states in sQGP. 

Any model of J /if) suppression must account for the 
experimental observation that a substantial fraction of 
the measured J/t/i's are from decay of the excited char- 
monium states Xc and -i/;/ [2l| . To calculate the survival 
probability of the excited states Xc and ■0', we use the 
same procedure as described above for the ground state 
J/V'. Above T^^ and T^,,, the excited states Xc and ij^' are 
assumed to melt. Further noting that the lattice moti- 
vated calculation indicate Tv 



~ T^i , we define a common 
above which all the states 
Xc and i/j' are dissolved. For feed-down fraction F, the 
J/tp survival probability is then obtained as, 



temperature — T^^ — T^,, 



^dis{T) = oo; T > Tj/^ 

Td^siT) = 0; T<Tj/^ (3) 

In EqlSl decay width of J/ip abruptly changes from 
oo to at the melting temperature Tj/^. It neglects 
the broadening of charmonium states below the critical 
temperature Tc- In a hot pion gas widths of decay and 
dissociation channels of the charmonium states ■0/, Xc 
and J/ip to DD, D*D, DD*, D*D* pairs are enhanced 
[2^. Lattice QCD calculations 24j also indicate D- and 



Sqgp = (1 - F)Sj/^ + FS^ (4) 

As mentioned earlier, PHENIX collaboration in d+Au 
collisions has studied cold nuclear matter effect on J/tp 
suppression [1,01 • J/i^'s are suppressed in d+Au colli- 
sions also. The suppression is consistent with Glauber 
model of nuclear absorption with J/-0-nucleon absorp- 
tion cross-section aabs = 2 ± 1 mb. If cold nuclear matter 
effect is taken into account, the survival probability of 
J/ip can be obtained as, 
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Sj/-4> — Sqgp X ScNAi, (5) 

where Scnm is the survival probabiUty in cold nuclear 
matter calculated in a Glauber model. 



III. RESULTS 

A. J/tp suppression in Au+Au collisions 

PHENIX data [1, on the centrality dependence of 
J/f/; suppression are shown in FiglS) In all the central- 
ity ranges of collisions, J/tp's are suppressed, suppres- 
sion increasing with collision centrality. One also notes 
that around Npart ~150, there is a distinct change of 
slope in the suppression. The black solid line in Fig[3] 
is an estimate of cold nuclear matter effect, calculated 
in the Glauber model with CTabs=2 mb. Experimen- 
tally, J/'i/i's are more suppressed than in the Glauber 
model calculation. Evidently, data demand suppression 
in addition to the nuclear absorption. In FiglH cen- 
trality dependence of J/ip suppression in the present 
"Hydro-|- J/f/;" model, for a choice of melting tempera- 
tures Tj/^/rc=l. 2, 1.4,1.6,1.8,2.0 and 2.2, are shown. As 
expected, for Tj/^ = 2.2Tc, J/tp are not suppressed in 
the QGP medium, the fluid temperature is below the 
melting temperature. J/tp's are increasingly suppressed 
as the melting temperature is lowered. It is apparent 
from Fig[31 that melting of J/ijj alone cannot explain the 
data. 

The "Hydro-|- J/?/'" model has three parameters, the 
melting temperatures, Tj/^ and and the fraction F 
of higher states {xc+i^')- If include the cold nuclear 
matter effect, then JZ-^-nucleon absorption cross-section 
aabs, can be considered as an additional parameter. How- 
ever, as mentioned earlier, cold nuclear matter effect in 
d-|-Au collisions is consistent with aabs = 2±1 mb. In the 
following, we assume that in Au+Au collisions also, cold 
nuclear matter effect is also adequately explained in the 
Glauber model of nuclear absorption with (Tabs =2 mb. 

Since the number of PHENIX data points are few, we 
do not fit simultaneously all the three parameters, Tj/^, 
Ty_ and F. For fixed melting temperature Ty,=l.-1.2Tc, 
we vary the Tj/^ and F to fit the PHENIX data. Note 
that in the present work Tc=164: MeV. The best-fitted 
values are shown in table 1. Melting temperatures Tj/^ 
and are well determined, Tj/^/Tc = 2.08 ± 0.25, 
T^/Tc = 1.1 ± 0.1. The feed-down fraction however 
can only be determined with large uncertainty, F = 
0.32 ± 0.20. Quality of the PHENIX data is far from sat- 
isfactory, data points are few, and the associated error 
bars are large. If quality of the data is improved, the un- 
certainty in F (and also in Tj/^ and T^) can be reduced. 
The melting temperature Tj/^ and and also the feed- 
down fraction F are in agreement with the results ob- 
tained by Gunji et o/.[l^. In FigHl the fit obtained to 



the data for T^=l.lTc, Tj/^=2.08Tc and F = 0.32 is 
shown. Sequential melting of charmonium states Xci "0' 
and J /ip in an expanding QGP medium, well explain the 
PHENIX data. 

Sequential melting of charmonium states in a decon- 
fined medium can explain the PHENIX data even if we 
disregard the cold nuclear matter effect. In table HI we 
have tabulated the best-fitted melting temperature Tj/^ 
and the feed-down fraction F obtained by fitting the data 
when cold nuclear matter effect is neglected. (crahs=0). 
Whether the cold nuclear matter effect is included or not, 
within the error, the melting temperature Tj/^ remain 
same. However, the data require higher feed-down frac- 
tion F « 0.5 for equivalent fit. In FigUl we have shown 
the fit obtained to the data with Tj/^ = 2Tc, = l.lTc 
and F — 0.5. The fit is comparable to the one obtained 
with the cold nuclear matter effect included. The re- 
sult is not unexpected. In ^1^ it was shown that in the 
phenomenological threshold model, sequential melting of 
charmonium states can explain the PHENIX data, with 
or without the cold nuclear matter effect. When cold 
nuclear matter effect is neglected for, the data demand 
higher feed-down fraction. Cold nuclear matter effects 
can be effectively mimicked by melting of higher states 
XcS and ijj' . The result underlies the importance of J/^p 
measurements in pA collisions. In heavy ion collisions, 
the feed-down fraction F can be accurately estimated 
only when the cold nuclear matter effect is accurately 
quantified in pA collisions. 

The present analysis suggests that the centrality de- 
pendence of J/^p suppression in mid-rapidity Au-f Au col- 
lisions at RHIC is consistent with sequential melting of 
Xc, '0' Eind J ftp in an expanding QGP medium. The 
melting temperature are well constrained, T^^ = T^, = 
1 - l.lTc, Tj/^ = 2- 2.1Tc. The feed-down fraction 
depends on whether the cold nuclear matter effect is in- 
cluded or not. Model require F — 0.3 ±0.2 if cold nuclear 
matter effect is quantified in the Glauber model with 
J/?/;-nucleon absorption cross-section aabs='^ mb. Larger 
fraction F = 0.5 ± 0.20 is required if cold nuclear matter 
effect is neglected. 

In all the calculations presented here, the decay width 
for J/ip abruptly changes from cx3 to at the melting 
temperature Tj/^ and similarly for the states Xc and ipf. 
Gunji et aZ.ll§| studied the effect of smoothening the de- 
cay width by using, 



^dis{T) oo; T > Tj/^ 

Td^siT) = a(T/re-l)2; T<Tj/^ (6) 

where a is the thermal width of the state at T/Tc=2. 
NLO perturbative calculations suggest that a > OAGeV 
[Hi. However, for a > 0.4, the PHENIX data are not well 
described. Data require that a < 0.1, when smoothening 
effect is not large. 
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FIG. 5: (color online) The filled circles are PHENIX data 
on the centrality dependence of J/tj} suppression in Cu+Cu 
collisions. The blue line is the absorption in Glauber model 
of nuclear absorption with aabs=2 mb. The solid line is the 
suppression due to sequential melting of Xc, ip' and J/ip in an 
expanding QGP fiuid, including the nuclear absorption. The 
dashed line is the same without the nuclear absorption. 



B. 



J/tp suppression in Cu+Cu collisions 



Recently PHENIX collaboration published their anal- 
ysis of J/ip measurements in Cu+Cu collisions [ll[. In 
Figini PHENIX data [ll| on the centrality dependence 
of J/tp are shown. In the most central collision, J ftp's 
are suppressed by a factor of ~ 2. Survival probabil- 
ity of J ftp increases as the collision centrality decreases 
till iVpart=16.7. In more peripheral {Npart=G-4:), colli- 
sions, survival probability decreases again. The decrease 
in survival probability as participant number decreases 
from to 6.4 is interesting. All the theoretical 

models predict continuous increase of suppression as the 
collision centrality increases. It was also noted [ll| that 
at comparable participant number, J ftp are suppressed 
similarly in Au+Au and in Cu+Cu collisions. However, 
there is a major difference between J ftp suppression in 
Cu+Cu and in Au+Au collisions. In FigO the blue line 
is the suppression in Cu+Cu collisions, calculated in a 
Glauber model of nuclear absorption, with (Tats =2 mb. 
As noted by the PHENIX collaboration, suppression in 
peripheral collisions (excluding Npart=&-^) is consistent 
with the nuclear absorption alone. In contrast, J/ip sup- 
pression in peripheral Au+Au collisions are not consis- 
tent with Glauber model of nuclear absorption. 

For hydrodynamical evolution of QGP fluid in Cu+Cu 
collisions, we use EqUto calculate the initial energy den- 
sity in the transverse plane. The constant £o in EqUJ 
depend only on the collision energy, not on centrality of 
the collisions. At the initial time Ti=0.6 fm/c, the initial 



FIG. 6: Initial temperature T{x,y — 0), in unit of the crit- 
ical temperature, of the QGP fiuid in Cu+Cu collisions at 
impact parameter b= 8.83, 7.55,6.82,6.00,5.05,3.87 and 2.29 
fm (bottom to top), respectively. 



energy density is e «30 GeV/fm^ or initial central en- 
tropy density is Sini — llOfm~^. In Fig. El the initial 
temperatures T{x,y — 0) of the fluid in b=2.29, 3.87, 
5.05, 6.0, 6.82, 7.55 and 8.83 fm Cu+Cu collisions are 
shown. The impact parameters roughly corresponds to 
participant numbers -/Vpart=98.2, 73.6, 53.0, 37.3, 25.4, 
16.7 and 6.4 respectively. In all these collisions, peak 
temperature of the fluid exceeds the melting temperature 
= T^^ — T,p,=l — l.lTc. Consequently, in all the cen- 
trality ranges of collisions, melting of the states Xc and ip' 
will contribute to the J/ip suppression. Immediately we 
find that J/4> suppression in peripheral Cu+Cu collisions 
is not consistent with sequential melting of charmonium 
states. Peripheral collisions are consistent with cold nu- 
clear matter effect and cannot accommodate additional 
suppression due to higher states x and Tpf. 

It is clearly shown in FigO When cold nuclear matter 
effect is included, PHENIX Au+Au data are explained 



(Fig|4]) with melting temperatures , T 



= l.lTc and feed-down fraction F 



J/i> 



2.08Te, 



0.3. In FigEl 

the solid line is the centrality dependence of J/ip sup- 
pression in Cu+Cu collisions, with cold nuclear matter 
effect included. The parameter values are unchanged, i.e. 
Tj/^ ^ 2.08Tc, = l.lTc and F = 0.3. As expected, 
the Hydro+ J/?/; model produces more suppression than 
required by the data. 

Will the data will be explained without any nuclear 
absorption but with increased feed-down fraction? As 
shown earlier, the PHENIX Au+Au data on J/tp sup- 
pression are explained with same melting temperatures 
but with increased feed-down fraction, F « 0.5. In FiglSl 
the dashed line is the centrality dependence of J/tp sup- 
pression without any nuclear absorption, but with in- 
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creased feed-down fraction F = 0.5. The data are not 
explained either. In contrast to Au+Au coUisions, se- 
quential melting of charmonium states in an expanding 
QGP medium do not explain the centrality dependence 
of J/ip suppression in Cu+Cu collisions. It may be noted 
that the present model neglect the recombination of cc 
pairs. Inclusion of recombination effect will reduce J/^ 
suppression. A shown earlier, within the present model, 
J/4' suppression in Au-t-Au collisions do not require any 
recombination. It is unlikely that recombination effect 
will be important in Cu-|-Cu collisions but not in Au-|-Au 
collisions. It appears that, mechanism of J/ip suppres- 
sion in Au+Au and in Cu-|-Cu collisions is different. 




-] — — — — 

1 10 100 1000 10000 

cm. energy (GeV) 

FIG. 7: Filled circles are the PHENIX data for the charged 
particle multiplicity per participant pair -^j^ — - ^ as a func- 
tion of cm. energy for participant number A'j,art=350. The 
solid line is a fit to the PHENIX data by EqU The unfilled 
circle is the extrapolated value of — at LHC energy 

-^3=5. 5 TeV, for participant number Npart = 350. 



IV. 3 111, SUPPRESSION IN PB+PB COLLISION 
AT LHC 

At the Large Hadron Collider (LHC), it is planned to 
collide Pb ions at centre of mass energy -v/s=5500GeV. 
Signals of the deconfinement phase transition are ex- 
pected to be better defined at LHC than at RHIC col- 
lisions. One wonders whether in LHC energy collisions, 
J ji^j will be more suppressed than at RHIC energy? Re- 
cently, in the phenomenological threshold model, it was 
predicted that suppression pattern in Pb-|-Pb collisions 
at LHC energy will be similar to that in Au-|-Au colli- 
sions at RHIC The reason is understood. In the 
threshold model, suppression depends on the local trans- 
verse density. If the local transverse density exceed a 



critical value, J/i/^'s are suppressed. The local transverse 
density depend marginally on the collisions energy and 
remain essentially same at RHIC and LHC energy col- 
lisions. In the "Hydro+ J/V'" model, J/i/' suppression 
depend directly on the "local fluid temperature" and in 
LHC energy collisions, local temperature will be substan- 
tially greater than that in RHIC energy collisions. One 
expects enhanced suppression at LHC energy. 




0.0 ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' 
50 100 150 200 250 300 350 400 

Npart 

FIG. 8: Hydro-|- J/i/) model predictions for the centrality de- 
pendence of J/V' suppression in Pb-|-Pb collisions at LHC en- 
ergy. The melting temperature for the states J/V', V'' ^-^d Xc 
are Tj/^ = 2.\T^ and = T^, = T^, = l.lTc The solid line 
is obtained with feed-down fraction F = 0.3 and including the 
nuclear absorption eff'ect with dabs =2 mb. The dashed line is 
without the nuclear absorption effect but with increased feed- 
down fraction F=0.5. For comparison, PHENIX data on the 
centrahty dependence of J/t/; suppression in Au+Au collisions 
are also shown. 

For predicting J/^/; suppression in Pb+Pb collisions at 
LHC, we have to guess the initial energy density distri- 
bution of the QGP fluid. PHENIX collaboration has 
tabulated the average charged particle multiplicity as a 
function of collision energy for a range of collision cen- 
trality. In Fig[71 for participant number A'pQj.t=350, the 
average multiplicity — T^^^^ shown as a function 
of collision energy . The multiplicity increases logarith- 
mically with energy, 

dMrh ^ ^ , I— , ^ 

—^^A + B\n^fs, 7 
drj 

with A = -0.33 and B = 0.75. We use the relation 
to extrapolate to LHC energy -y/s=5.5 TeV. The extrap- 
olated value of average charged particle multiplicity in 
LHC energy is ~ 927 ± 70. We adjust the central en- 
tropy density to 5'i„i=180 fm~^ such that a Npart=-^^0 
Pb+Pb collisions produce ~ 900 charged particles. The 
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participant numbers in Pb+Pb collisions are calculated 
with NN inelastic cross-section crinei=70 mb. Entropy 
density Sini = 180 fm~^ corresponds to central tempera- 
ture corresponds to Ti— 421 MeV. This can be contrasted 
to central temperature Ti—357 MeV in Au+Au collisions 
. Compared to Au+Au collisions at RHIC, in Pb+Pb 
collisions at LHC, central temperature approximately ~ 
20% higher. 



3.0 



(a) Au+Au 
Tc=164 MeV 



(b) Pb+Pb 
Tc=164 MeV 




-8 -6 -4-2 2 4 6 
X(fm) 



6 -4 -2 2 4 6 
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FIG. 9: In panel (a) and (b), initial temperature T{x,y = 0) 
of the QGP fluid, in unit of the critical temperature (Tc), in 
Au+Au collisions at RHIC and in Pb+Pb collisions at LHC, 
are compared. The lines (top to bottom) corresponds to 0- 
10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60-70%, 70- 
80% collisions. 

In FiglH Hydro+ J/-!/) model predictions for the central- 
ity dependence of J /ip suppression in Pb+Pb collisions 
are shown. The melting temperature of the states J /ip 
and Xc and ipt are property of the deconfined medium 
and do not depend on the collision energy. We use 
the values extracted from the analysis of RHIC data. 



70%, 70-80% and 80-90% in Au+Au colhsions RHIC and 
in Pb+Pb collisions at LHC. Both in Au+Au and Pb+Pb 
collisions, peak temperature in 0-10% centrality collisions 
exceed the melting temperature Tj/^ = 2.1Tc. In 0- 
10% centrality collisions, J/ip will be suppressed both 
in RHIC and LHC energy collisions. However, it is not 
so in less central collisions, e.g. in 20-30% centrality col- 
lisions, while peak temperature exceed the melting tem- 
perature Tj/^ in Pb+Pb collisions at LHC, in Au+Au 
collisions, the peak temperature is less than the melt- 
ing temperature. J/%1) will survive in 20-30% centrality 
Au+Au collision at RHIC but they will be melted in 20- 
30% centrality Pb+Pb collisions at LHC. Thus J/ip's can 
survive in mid-central Au+Au collisions at RHIC, but 
are suppressed in mid-central Pb+Pb collisions. For the 
very reason, the change of slope in the suppression pat- 
tern also occurs at lower participant number. Compared 
to RHIC energy collisions, the melting temperature Tj/^ 
is reached in lower participant number collisions at LHC 
energy. 



V. SUMMARY AND CONCLUSIONS 

To summarise, we have developed a "Hydro+ J/?/;" 
model to study J/^p suppression in an expanding Quark- 
Gluon-Plasma. The space-time evolution of the QGP 
fluid is obtained by solving the hydrodynamic equations 
for ideal fluid in 2+1 dimensions. At the initial time, 
J/ij/s are randomly distributed in the transverse plane. 
As the fluid evolve in time, the free streaming J/ip's are 
completely suppressed if the local fluid temperature ex- 
ceed a critical temperature Tj/^. Similarly, the states 
Xc and ipi are assumed to melt above a critical tempera- 
ture T^. The melting temperatures Tj/^ and and the 
feed-down fraction F from the higher states Xc and ■0/, 
are fitted to reproduce the PHENIX data on the central- 
ity dependence of J /ip suppression in Au+Au collisions 
at mid-rapidity. The PHENIX data are well explained 



T 



2 - 2.1Tc and = I - I.IT^. In Figli the with Tj/.^ « 2 - 2.1Tc, ^ I - I.IT^. The feed-down 



solid line is the suppression obtained with Tj/^ = 2.irc 
and = l.lTc, F — 0.3, the cold nuclear matter ef- 
fect is included. Nearly similar suppression is obtained 
when cold nuclear matter effect is neglected but, the feed- 
down fraction is increased to F = 0.5 (the dashed line in 
FiglH]). For comparison, we have shown the PHENIX 
data on J /tp suppression in Au+Au collisions. In very 
central collisions or very peripheral collisions, J /ip are 
suppressed similarly in RHIC and LHC energy. But in 
mid central collisions, J/4>^s are more suppressed in LHC 
energy than in RHIC energy collisions. One also notices 
that while in Au+Au collisions at RHIC, the change in 



slope occur around Nj 



part 



150, in Pb+Pb collisions 



at LHC, the change occur at lower participant number 
Npart ~ 70. The reason can be understood from Fig. [51 
where we have compared initial temperatures T(x, y = 0) 
in 0-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60- 



fraction depends on whether or not J/^p suppression in 
cold nuclear matter is included. If cold nuclear matter 
effect is included, data require F w 0.3. Fraction F in- 
creases to F wO.5, if the cold nuclear matter effect is 
neglected. It appears that, to a large extent, suppression 
of the states Xc and ^f, can mimic the cold nuclear matter 
effect. While sequential melting of Xc and J/ip in an 
expanding QGP fluid well explain the J/ip suppression 
in Au+Au collisions, the model fails to reproduce the 
experimental data in Cu+Cu collisions. In Cu+Cu col- 
lisions, "Hydro+ J/t/i" model produces more suppression 
than required by the data, indicating that J/ip suppres- 
sion mechanism in Au+Au and in Cu+Cu collisions are 
different. We have also given prediction for the central- 
ity dependence of J/tp suppression in Pb+Pb collisions 
at LHC energy. The model predicts more suppression in 
LHC than in RHIC energy collisions. 
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TABLE I: Melting temperatures (in unit of Tc) of direct J/ij}, 
{Tj/nj,), the states Xc s^nd V'' (^Ix:) and the fraction of the 
higher states Xc and tp' . For fixed melting temperature T-^, the 
PHENIX data are fitted by varying the melting temperature 
Tj/^ and F, Values obtained with (cra63=2 mb) and without 
{o'abs= mb) the cold nuclear matter effects are shown. 



Tr 


crabs='2 mb 

F xVd.o.f 


O"ab3=0 

^ F xVF 


f.OO 


2.08 ± 0.18 0.30 ±0.17 0.72 


2.07 ±0.22 0.49 ±0.22 1.50 


1.05 


2.08 ± 0.23 0.31 ±0.18 0.72 


2.00 ±0.22 0.51 ±0.18 0.89 


1.10 


2.08 ± 0.13 0.32 ±0.18 0.68 


2.00 ±0.33 0.51 ±0.18 0.90 


1.15 


2.07 ± 0.46 0.32 ±0.23 0.97 


2.05 ±0.37 0.55 ±0.25 1.48 


1.20 


2.08 ± 0.27 0.32 ±0.23 0.96 


2.00 ±0.51 0.56 ±0.26 1.47 
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